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ABSTRACT 

We are conducting a survey for pulsars and transients using the Giant Me- 
trewave Radio Telescope (GMRT). The GMRT High Resolution Southern Sky 
(GHRSS) survey is an off-Galactic-plane (|6| > 5) survey in the declination range 
—40° to —54° at 322 MHz. With the high time (up to 30.72 /rs) and frequency 
(up to 0.016275 MHz) resolution observing modes, the 5a detection limit is 0.5 
mJy for a 2 ms pulsar with 10% duty cycle at 322 MHz. Total GHRSS sky cov¬ 
erage of 2866 deg^, will result from 1953 pointings, each covering 1.8 deg^. The 
lOn detection limit for a 5 ms transient burst is 1.6 Jy for the GHRSS survey. 
In addition, the GHRSS survey can reveal transient events like rotating radio 
transients or fast radio bursts. With 35% of the survey completed (i.e. 1000 
deg^), we report the discovery of 10 pulsars, one of which is a millisecond pul¬ 
sar (MSP), one of the highest pulsar per square degree discovery rates for any 
off-Galactic plane survey. We re-detected 23 known in-beam pulsars. Utilising 
the imaging capability of the GMRT we also localised 4 of the GHRSS pulsars 
(including the MSP) in the gated image plane within ± 10". We demonstrated 
rapid convergence in pulsar timing with a more precise position than is possible 
with single dish discoveries. We also exhibited that we can localise the bright¬ 
est transient sources with simultaneously obtained lower time resolution imaging 
data, demonstrating a technique that may have application in the SKA. 
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1. Introduction 


Pulsars are rapidly rotating neutron stars emitting a beam of radio waves from their 
magnetic poles, that sweep our line-of-sight at the spin period. In the 48 years since the 
discovery of pulsars, the present population of about 2300 known pulsars have spin periods 
ranging from 1.4 ms to 8.5 s, magnetic held strengths from 6.6x10^ G to 1.2x10^^ G, 
line-of-sight electron column density values, called dispersion measures (DMs), from 2.38 
pc cm“^ to 1456 pc cm“^, hux densities at 400 MHz ranging from 0.1 mJy to 5000 mJy 
and characteristic ages ranging from 2.2x10^ years to 6.7x10^° year^ In an attempt to 
understand the distribution of sources within the limits of the known population and to 
probe the limits themselves there are many ongoing surveys that are discovering pulsars at 
an encouraging rate. Although many of the newly discovered pulsars have similar properties 
to the known population, some of the new ones are pushing the boundary of known parameter 
space. For example Low-mass X-ray binary—radio millisecond pulsar transitioning system: 


J1023-I-0038 (Archibald et al. 2009 

) and J1227-4853 

Roy et al.|2015), a hierarchical stellar 

triple pulsar system: J0337-I-1715 ( 

Ransom et al.||2014 

1 , a massive neutron star—white dwarf 

binary: J0348-I-0432 (Antoniadis et al.||2013; 

, and a pr 

dsar with a very dense planetary mass 
iversity of properties of these objects 

companion: J1719—1438 (Bailes et al. 2011 

). The d 


justihes the requirement of having more sensitive surveys for further discoveries. 


Modeling by Faucher-Giguere et al. (2006) indicates that the Galaxy contains about 10° 
pulsars and the number of presently known pulsars is only around 1% of this population. 
This implies that a vast majority of pulsars are waiting to be discovered. Studies of pulsars 
yield a better understanding of a variety of physics problems, from acceleration of particles 
in ultra strong magnetic helds (primarily via study of emission properties of normal pulsars, 
having spin period > 30 ms) to probes of ultra dense matter (mostly via studying the timing 
properties of millisecond pulsars, having spin period <30 ms, that are very stable rotators). 
Pulsars provide useful probes of their environments, e.g. inside pulsar wind nebulae, the 
centre of Globular clusters, or the Galactic centre. Some of the normal pulsars are useful for 
investigation of single pulse behaviour and can exhibit interesting individual properties, like 
glitches (Lyne et al. 1996), prohle state changes (Lyne et al. 2010), nulling (Backer 1970) 


and intermittency (Kramer et al. 2006). The emission mechanism for pulsar radiation at 


radio wavelengths is not yet understood, and a range of theoretical models have been pro¬ 
posed (e.g. Ruderman et al. (1975); Gil et al. (2003); Alexander et al.| (2014)) to explain it. 
The discovery of a new sample of pulsars with diverse observational properties will provide 
additional constraints on the possible emission mechanism. Additional pulsars distributed 


® ht tp: / / WWW. atnf. csiro. au/people/pulsar/psrcat / 
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in our Galaxy will aid the investigation of properties of the interstellar medium via scat¬ 
tering measurements as well as via dispersion measure and rotation measure studies. The 
fractional rotational stability of millisecond pulsars (MSPs), one part in 10^^, is comparable 
to atomic clocks ( Lorimer et al.||2004 ). Such rotational stability, compactness second only to 
black holes, and their presence in binary systems, make MSPs ideal laboratories to test the 


physics of gravity and as detectors for long-wavelength gravitational waves (Detweiler et ah 


1979). Some relativistic binary MSPs are useful for tests of gravity and an array of spatially 


distributed MSPs can be used to detect gravitational waves (Lee et ah 2012). In addition 


to the contribution towards the detection of gravitational wave signal, MSP evolutionary 
processes can be tracked through individual interesting discoveries of MSPs in special evo¬ 


lutionary phases. For example, the black-widow systems (Roberts et ah 2011) provides a 


missing link between the binary and isolated MSPs and Low-mass X-ray binary—radio MSP 


transitioning redback systems (Archibald et ah 2009 Roy et ah 2015) provides a probe of 
binary evolution. 

In addition to the regular emission from pulsars, a wide variety of transient phenomena 
are known at faster time scales (micro-seconds to seconds). The discoveries of giant pulses 


from pulsars (Lundgren et ah 1995), quasi-periodic emission from rotating radio transients 


(RRATs; McLaughlin et al. (2006)), bursts of periodic pulsations from magnetars (Camilo 


et al. 


and highly dispersed fast radio bursts (FRBs) of possible extragalactic origin 


(Lorimer et al. 2007) are a few examples. The FRBs have been discovered in radio pulsar 


surveys at Parkes and Arecibo in the last decade ( 

Lorimer et al. 

2007 

Thornton et al. 

2013 

Spitler et al. 

2014, 

Burke-Spolaor & Bannister et ah] 

2014 

Petroff et al. 

2014 

). All FRBs 


discovered to date have been single radio events of millisecond duration with DM values 
generally higher than the possible Galactic contribution. Majority of the FRB population 
theories suggest that they originate at cosmological distances (Deng et al. (2014|); Luan et 


al. (2014); Keane & Petroff (2015)), although the actual progenitors are still unknown. The 


probable extragalactic origin allows these bursts to be used for determination of the baryon 
content of the intergalactic medium. The detection of a large number of FRBs (presently 
only 17 FRB^ are known) is required in order for us to have a better understanding of the 
nature of the population. This will allow one to establish whether they form a population 
of standard candles or some other distribution. Also many lines-of-sights will be needed to 
make them useful probes of the intergalactic medium and for determining the presence of 
the so-called missing baryons. With a larger population of sources it will also be possible to 
determine the spectral index of the radio emission which is a key element for probing the 


emission physics. Magnetar flares (Thornton et al. (2013); Kulkarni et al. (2014)), blitzars 


^https://astro.uni-bonn.de/^tauris/NS2015/Keane_FRBs.pdf 
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(Falcke et al. 2014), super-giant pulses from neutron stars (Cordes et al. 2015), and pulsar- 


planet systems (Mottez & Zarka 2014), are some of the models suggesting extragalactic 


origin of FRBs at cosmological distances. 

The known pulsar population is increasing with the surveys running at major single 
dish and array telescopes around the world: e.g. High Time Resolution Universe (HTRU) 


survey at Parkes (Keith et al. 2010), HTRU-North survey at Effelsberg (Barr et al. 2011), 


survey for Pulsars and Extragalactic Radio Bursts (SUPERB) surveyn at Parkes, Pulsar 


survey Arecibo L-band Feed Array (PALFA, Cordes et al. (2006)), Green Bank Telescope 


drift scan survey (Boyles et al.||2013). Green Bank Northern Celestial Cap (GBNCC) survey 


(Stovall et al. 2014), Arecibo all-sky 327 MHz drift pulsar survey (A0327, Deneva et al. 


(2013)), the LOFAR Pulsar Pilot Survey (LPPS) survey (Coenen et al. 2014), the LOFAR 


Tied-Array All-Sky Survey (LOTAAS)|^ at LOFAR and Fermi directed MSP surveys (Ray 


et al. 2012). 


Because of the generally steep spectral nature of pulsars, lower frequencies are an obvious 
choice for searching for fainter pulsars away from the Galactic plane, where search sensitivity 
is not severely affected by sky temparature and increased scattering. Such surveys away from 
the Galactic plane will detect relatively older pulsars. The success of Fermi directed searches 
in discovering more than 60% of Fermi MSPf at frequencies below 1 GHz demonstrates the 


advantages of low frequency searches ( |Ray et ^|2012 ). Besides the GBT and the LOFAR, 
the Giant Metrewave Radio Telescope (GMRTp^ is another facility which has a sensitive 
low frequency (< 600 MHz) observing capability. With the aid of reduced quantised noise 


supported by the flexible GMRT Software Backend, (GSB; Roy et al. (2010)) the search 
sensitivity to pulsars in general is improved by 30% compared to its old hardware counterpart. 
An additional 30% increase in sensitivity to MSPs is achieved through using the higher 
time resolution mode implemented in the GSB. Thus, effectively, the search sensitivity to 
normal pulsars is improved by 30% whereas for MSPs a 60% improved sensitivity is achieved 


with the GSB ( 

Roy et al. 

2013 

). This is exemplihed by the discovery of seven MSPs from 

Fermi directed searches by 

Bhattacharyya et al. 

(2013 

). These were the hrst Galactic MSPs 


discovered with the GMRT. This clearly illustrates the potential of low frequency pulsar 
search using the GMRT. 


®https://sites.google.com/site/publicsuperb/ 

® ht t p: / / WWW. ast ron. nl/ lot aas 

^°http: / / astro.phys.wvu.edu/GalacticMSPs / GalacticMSPs.txt 
^^http://gmrt.ncra.tifr.res.in 
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Since we have no information on most pulsar positions a priori (unlike the ones associated 
with unidentihed Fermi sources or supernovae remnants), blind searches are required for 
pulsar discoveries. Lower frequency surveys beneht from a larger £eld-of-view, in addition to 
the steep spectral nature of pulsars. The success of targeted searches with the GMRT makes 
it evident that blind searches with the GMRT are warranted and can produce a signihcant 
science yield in terms of discovering new pulsars, MSPs and transient events. Thus we 
are carrying out the GMRT High Resolution Southern Sky (GHRS^^ survey for pulsars 
and transients to conduct one of the most sensitive and highest resolution surveys for the 
southern hemisphere. To make the search more efficient and complementary we target a 
part of the southern sky which has not been searched at frequencies below 1.4 GHz for the 
last two decades. New observing modes were developed to increase the time and frequency 
resolution of the survey, making it even more sensitive to MSPs and normal pulsars with 
higher DM values. 

Large uncertainties associated with the discovery positions (e.g. ± 40' for GMRT—322, 
± 18' for GET—350, ± 7' for Parkes—Lband) hinder sensitive follow up studies of these newly 
discovered pulsars using coherent beams of array telescopes, or at higher frequencies using 
single dishes. The GMRT interferometric array allows us to localise the newly discovered 
pulsars and transients in the image plane with an accuracy of better than ± 10" (half of the 
typical synthesized beam used in the image made at 322 MHz). Precise a priori astrometric 
positions are also needed to overcome the effect of large covariances in the timing hts (with 
discovery position, pulsar period derivative and unknown binary model in case of binary 
pulsars) for rapid convergence of an initial timing model. Moreover, sensitive coherent array 
follow up observations can significantly reduce the use of telescope time by ~ 20x for the 
GMRT (due to the 4 to 5 times sensitivity improvement for coherent array observations) and 
also improves the uncertainties in time-of-arrivals (TOAs) or generates more closely spaced 
TOAs in order to avoid ambiguous phase connection while timing. 

Finding pulsars is only half the job. The other half is to understand their nature, 
rotation properties, spectral indices and relation to the pulsar population; for this regular 
follow up is needed. Precise measurements of these parameters puts a new discovery in its 
proper place in the pulsar population. Equipped with the precise localisation we carry out 
more sensitive follow up observations with the coherent array of the GMRT. In this paper 
we demonstrate that convergence in the timing £t is facilitated by a priori knowledge of a 
precise position. 

All the known FRBs detected to date have been found with single dish telescopes, and 


^^http://www .jb.man.ac.uk/research/pulsar/Resources/ghrss.html 
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for most of them the uncertainty in position is of order of ± 7' or more. Interferometric 
imaging using the simultaneously recorded visibility data at 2 s interval will allow precise 
localisation of the FRBs detected from the GHRSS survey. This will allow multi-wavelength 
follow up and thus the possibility of identification of putative host galaxies. However, this 


will only be feasible for the relatively stronger (> 3 Jy) (detailed in Section 9.2) FRBs. Thus 


we will be able to localise bright FRBs (e.g. Lorimer bursts, Lorimer et ah (2007)). 


In this paper we describe the survey, the data analysis pipeline, discovery parameters of 
the 10 pulsars, the localisation of 4 and timing solution of 2 pulsars discovered in the GHRSS 
survey. The observing system and target sky are detailed in Section and Section of the 
paper. The prediction of survey outcome is detailed in Section The survey processing is 
discussed in SectionIn Section]^ we calculate the theoretical sensitivity and compare that 
with the sensitivity achieved in the survey. Section describes the re-detection of known 
pulsars from this survey and their interferometric flux measurements. In Section we give 
the basic parameters of the pulsars discovered in the GHRSS survey and in Section we 
detail the localisation of detected events using the GMRT interferometric array. The timing 
and follow up of the pulsars discovered in the GHRSS survey is described in Section The 
relevance of GHRSS survey results and methods with the future Square Kilometre Array 


(SKA) is discussed in Section 11 Section 12 presents the summary of this paper 




0 


Fig. 1.— (a) Sky coverage in Galactic coordinates of recent and ongoing pulsar surveys 
around 300 MHz; GBT drift scan (pink) at 350 MHz, GBNGG (purple) at 350 MHz and 
GHRSS (green) at 322 MHz. Red dots are the pulsars discovered with the GHRSS survey, 
(b) Sky coverage in RA and DEG. 
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2. Observing system 


The GMRT is a multi-element aperture synthesis telescope consisting of 30 antennas, 
each of 45 m diameter, spread over a region of 25 km diameter and operating at 5 different 


wave bands from 150 MHz to 1450 MHz (Swarup et ah 1997). The GMRT Software Backend, 


built using GOTS components, is a fully real-time backend that supports a FX correlator 
and a beamformer for an array of 32 dual polarized signals, Nyquist sampled at 33 or 
66 MHz (Roy et al. 2010). We developed high time and frequency resolution observing 


modes that are ported to the GMRT Software Backend, producing incoherent or coherent 
beam 2048x0.016275 MHz hlter-bank outputs sampled every 61.44 ps (hereafter HiResl) or 
1024x0.03255 MHz sampled every 30.72 /rs (hereafter HiRes2). This translates to a data rate 
of 32 MB/s for the high time-resolution data with 8 bits/sample after adding two polarized 
intensities. In parallel we recorded frequency averaged 512 channel visibilities for each of 
the 465 baselines at 2 s intervals resulting in a 2 MB/s data rate. The target scans were 
interleaved with scans on nearby calibrator sources every 1 hour to allow for phase calibration 
for imaging. The visibility data were then converted to FITs hies for further processing in 
13 The GMRT interferometric array provides the opportunity for combining 


AIPg^l/CAS 


the imaging mode with beamforming to widen the scope of studying pulsars (e.g. Roy et 


al. (2012); Bhattacharyya et al. (2013); Roy et al. (2013)). For the GHRSS survey we used 
the incoherent array beamformer output in order to maximise the £eld-of-view. In order to 
optimise the gain of the incoherent beam, the antenna based gain offsets were calibrated out 
prior to the addition of the signals. 

Gombining contributions from receiver temperature, ground temperature and spillover 
from the ground, the total system temperature at 322 MHz is 106 K on cold sk} The data 


were recorded to local disks and then transferred to the dedicated 64 TB GHRSS storage at 
the National Gentre for Radio Astrophysics (NGRA) using an 8 Gbps hbrelink. We shipped 
the data disks from the GHRSS storage to the University of Manchester for processing on 
the Hydrus supercomputer. The Hydrus supercomputer with a dedicated 64 TB storage was 
used as a major processing and data-centre host for the GHRSS survey. Raw data were 
archived on LT05 magnetic tapes immediately after the observations for long term storage. 


^^http://www.aips.nrao.edu/index.shtinl 
^"^http: / / casa.nrao.edu 

^^http://www.gmrt.ncra.tifr.res.in/gmrt_hpage/Users/doc/GMRT—specs.pdf 






















3. Sky coverage 


The 322 MHz GHRSS survey consists of scans of the sky (ranging from —40° to —54° in 
declination), complementary to other low frequency surveys at the GET, and LOFAR. Figure 
shows portion of the target sky (marked in green) compared to the GET surveys at 350 
MHz, noting that LOFAR is not scanning below declination zero. Though the pink and 
purple bands are indicative of GET drift scan survey and GENGG surveys, there are empty 
patches in the drift scan survey, which are being covered by the GENGG survey. 


The target sky for the GHRSS survey has not been searched at frequencies below 1.4 GHz 
for the last two decades. The Parkes southern sky survey at 436 MHz with ~3 mJy search 


sensitivity ( 

Eailes et al. 

1994 

was conducted from 1991 to about 1995, and it discovered 85 

normal pulsars and 17 MSPs 

16 



For the GHRSS survey we are not considering low Galactic latitudes (|6|<5°). The 
target sky has two components, a mid-latitude component, MGL, for 5<|6|<20° covering 
the sky with declination range —40° to —54° with 1200 s integration with HiResl mode, and 
a high latitude component, HGL, (|&|>20°) covering the same declination range with 900 s 
integration with HiRes2 mode. The range of sky temperature for the GHRSS survey range 
from 33 K to 220 K based on the Haslam et al. ( |1982| ) 408 MHz map. 


In order to partly compensate for the ~ 50% increase of sky temperature at MGL (|6|~ 
10°) with respect to HGL, we choose to increase the integration time to 1200 s. The target 
sky is covered with 682 and 911 pointings respectively for the MGL and HGL components, 
totaling 1953 pointing scans, each covering 1.8 deg^ adding up to a total coverage of 2866 
deg^. The survey parameters are summarised in Table 


4. Survey prediction 


In order to determine the discovery potential of the GHRSS survey we performed sim¬ 
ulations based on the PSRPOPP\f^ode ( |Eates et al. 2014). psrpoppy is a modihed version 
of PSRPOP (Lorimer et al.||2006[) that models the period, luminosity and spatial distributions 


of a set of pulsars discovered in the successful Parkes Multibeam Pulsar Survey (Manchester 


et al.||2001 ) and the high-latitude pulsar survey ( Eurgay et aL||2006 ) to predict the discovery 
potential of the ongoing surveys (e.g. Keith et al. (2010); Keane et al. (2015)). For the 


^®http://www.jb.man.ac.uk/pulsar/research/jodsum/node3.html^s70 
^^https: //github.com/sambSs/PsrPopPy 
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Fig. 2.— Integrated profile and phasogram of PSR B1007—47 for before (left panel) and 
after (right panel) zero-DM RFI mitigation. This observation was nnusually affected by RFI. 
Even with some residual RFI present, a signihcant cleaning of the data and improvement of 
SNR by about a factor of 3 is achieved. The zero-DM RFI mitigation is followed by further 
cleaning of the data using the rfifind package of presto resulting in 8%-10% improvement. 


purpose of simulation we considered normal pulsars defined as P > 30 ms, P < lx 10“^^. 
In simulation of the expected normal pulsar population from the GHRSS survey, we used 
the input parameters listed in Table (which are generally default options of psrpoppy and 
are similar to the ones used in Keith et ah (2010)). 

Using the parameters listed in Table [^and considering a total of 30,000 possibly detectable 
pulsars, the psprpoppy software generates a model ( Bates et ah [ 20 14) such that with Parkes 
Multibeam Pulsar Survey, 1038 pulsars are detected, which is close to the 1122 pulsars ac¬ 
tually discovered in this survey. With this model, we ran the simulation 100 times with 
the specifications of the GHRSS survey and predict that it would discover 97±7 pulsars, 
considering spectral index mean of —1.7. Simulations considering different spectral indices 
with mean of —1.4 and —1.8 result in discovery predictions of about 81±6 and 113±8 normal 
pulsars respectively. We expect to find considerably fewer MSPs due to their intrinsic faint 
nature in addition to dispersion smearing and scatter broadening. For the simulation of 
the MSP population we considered a modihed input for the spin period distribution, height 
above the Galactic plane and luminosity distribution as listed in Table According to the 
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Dispersion measure (pc cm 


Fig. 3.— Scatter broadening and dispersion smearing as a function of DM. At DM<75 
pc cm“^ scatter broadening dominates over dispersion smearing. Thus the GHRSS sur¬ 
vey sensitivity is not significantly affected by the instrumental broadening up to DM of 75 
pc cm“^. 


PSRPOPPY simulation the GHRSS survey will discover around 9±3 MSPs considering spec¬ 
tral index mean of —1.7. Simulations considering spectral indices with mean of —1.4 and 
—1.8 resulted in discovery predictions of about 7±3 and 11±4 MSPs respectively. 


We also tried to judge the discovery potential of the GHRSS survey by comparing with 
the ongoing GBNGG survey which is scanning the sky at a similar frequency and with 
comparable sensitivity to the GHRSS survey i.e. about 0.5 mJy of minimum detectable flux 
density (detailed in Section]^. GBNGG survey has covered 19500 deg^ and has discovered 
108 pulsar^^ Scaling from the GBNGG areal discovery rate, we expect to find of the order 
of 16 pulsars in the total GHRSS sky coverage (2866 deg^). 


In addition to the potential of discovering pulsars and millisecond pulsars, the GHRSS 
survey provides a snapshot of the transient radio sky in the southern hemisphere at 322 
MHz. Based on the 35% completion of HTRU-South survey analysis, Thornton et al. (2013) 
estimated the FRB rate as ij^QgXlO^ sky“^ day“^ for bursts having fluence > 3 Jy ms. For 


^®http://arcc.phys.utb.edu/gbncc/ 
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Fig. 4.— Top panel ; Example diagnostic plot from the GHRSS periodicity search pipeline. 
This plot shows the MSP J2144—5237 discovered in the GHRSS snrvey with incoherent array 
observations, Bottom panel; presto search ontpnt of PSR J2144—5237 with a coherent 
array observation. The plot demonstrates the improved detection signihcance with coherent 
array observations which can be performed after precise localisation of the pnlsar in the 
image plane (Section]^. 
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Single pulse results for ’HR_1559-44_!A_$earch.lodm40' 
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Fig. 5.— Example output of our presto based single-pulse search pipeline using CPU 
processing, showing bright single pulses from PSR J1559—4438 (which is a known pulsar 
re-detected in the GHRSS survey) around a DM of 56 pc cm“^. The top left panel shows the 
signal-to-noise distribution of the detected pulses, number of pulses versus trial DM (in top 
centre), and signal-to-noise as a function of trial DM. The lower plot shows signal-to-noise 
of events as a function of time and trial DM. 


the GHRSS survey we calculated a minimum detectable flux of 1.6 Jy for lOcx detection with 
5 ms bursts (detailed in Section]^. Goincidentally this minimum detectable flux for an FRB 
translates to a fluence of ~ 3 Jy ms. Gonsidering the total 2866 deg^ of the GHRSS survey 
area with ~ 15 mins of dwell time, we expect to discover FRBs having fluence > 3 Jy 
ms. However this prediction is based on the assumption that FRBs have a flat spectrum all 
the way from 1.4 GHz to 300 MHz and of course scattering and spectral index are largely 
unknowns for FRBs. Thus, the predictions are likely be lower than these estimates unless 
they are steep spectrum objects, and scattering at low frequencies is not large. Moreover the 
uncertainty on the rate presented by Thornton et ah (2013) and new lower rate estimates 
from Rane et ah (2015), Ghampion et al. (2015) and Petroff et ah (2014), the number of 
detected FRBs could be a factor of a few lower. However, any GHRSS detections will help 
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Fig. 6.— Example output of biforst single pulse search pipeline with GPU on the known 
pulsar J1559—4438 (re-detected in the GHRSS survey) showing bright pulses around a DM 
of 56 pc cm“^. The top left panel shows the detected pulses versus DM and top right panel 
shows signal-to-noise as a function of trial DM. The lower plot shows signal-to-noise of events 
as a function of time and trial DM. 


us to signihcantly improve rate determinations and the frequency dependence of FRBs. 


The comparison of survey predictions with achieved discovery rate of pulsars and tran¬ 
sients from the GHRSS survey is discussed in Section 


5. Survey processing 

To analyse the survey data, we used the 1456 core Hydrus cluster (30 Tflops) of the 
Jodrell Bank pulsar group and 512 core IBM cluster (10 Tflops) at the NGRA. The full 
GHRSS survey will generate approximately 50 TB of data and so far we have collected 
about 17 TB data corresponding to 35% completion of the survey. To optimise the data 
processing (i.e. considering the data transfer time from the observatory to the University of 
Manchester) and maximise the usage of the available GPU cycles on the available computing 
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clusters we analysed 50% of the data with Hydrus and the rest with the IBM cluster. 


5.1. RFI flagging 

Like most low frequency radio telescopes, the GMRT is affected by local RFI. Frequently 
observed RFI occurrences at the GMRT include broadband impulsive signals originated from 
spark discharges from power-lines, narrow spectral impurity and short time-domain bursts. 
We performed zero-DM subtraction (Eatough et ah 2009) to remove broad band interference. 
The zero-DM subtraction procedure signihcantly reduced the number of spurious candidates 
detected in periodicity searches and in single pulse searches. This dramatically improved 
the detection signiflcance towards broad band impulsive signals coming from pulsars. Fig¬ 
ure]^ is the integrated prohle and phasogram of PSR B1007—47 before and after zero-DM 
subtraction. Signihcant RFI can be seen as diagonal stripes caused by pulsed emission from 
power-lines or its subsequent harmonics, which are mitigated to a large extent after zero-DM 
subtraction and signihcant improvement of the pulsed detection is observed. Some residual 
RFI is still present in the data which can be seen in Figure rfifind generally cleans up 
the data, especially for near to zero-DM residual RFIs, narrow-band spectral features and 
short time-domain bursts. We made this zero-DM subtraction routine a part of the pub¬ 
licly available SIGPROC distribution. Fourier domain analysis of the time series dedispersed 
at a DM = 0 enabled us to excise the powerline frequency, 50 Hz. Finally, we employed 
PRESTCp^based rfifind package on the zero-DM subtracted time series to remove any re¬ 
maining bursts of interference or strong spectral features. It is noteworthy that without 
these RFI flagging steps, 90% of the discovered pulsars reported in this paper would not 
have been detected. 


5.2. Dedispersion 

We dedispersed the data over a range of DMs, 0 to 500 pc cm“^, resulting in about 6000 
dispersion trials for HGL and 10000 for MGL. The dispersion plan is generated following 
the DDPLAN of PRESTO in order to minimise the loss of sensitivity caused by dispersion 
smearing over a single frequency channel (Table |^. We reduced the time resolution for 
higher dispersion measure trials by down-sampling the data. For the GHRSS observing 
mode with 2048 frequency channels, the dispersive smearing within a frequency channel for 


^®http://www.cv.nrao.edu/~sransom/presto/ 





15 


a DM of 500 pc cm ^ is 1.9 ms at 322 MHz. Figure]^ shows the scatter broadening (|Bhat et 


al. [2004) and dispersion smearing for the GHRSS survey as a function of DM. This indicates 


that with the newly implemented high resolution mode, the GHRSS survey sensitivity is not 
signihcantly affected by the instrumental broadening resulting from dispersion smearing up 
to a DM of 75 pc cm“^. 


5.3. Periodicity search 

The dedispersed time series were Fourier transformed and candidates with periodicities 
of known RFI sources were removed. We used a Fourier-based acceleration search method us¬ 
ing the standard pulsar search techniques implemented in presto. The acceleration search 
spanned a broadening of the Fourier peaks by up to 50 Fourier frequency bins (i.e. Zmax of 
50). This translates to 5 m s“^ line-of-sight acceleration for a 2 ms pulsar over 15 mins of ob¬ 
serving scan detected in all 8 harmonics. The search for periodicity was done using harmonic 
summing (up to 8 harmonics). The candidates for each DM were stored for comparison. 


5.4. Candidate sorting 

Gandidates within a Fourier bin of another candidate were considered to be the same. 
We also ignored the candidates with summed incoherent power less than 4cr. Harmonically 
related candidates were removed, keeping the one with highest detection signihcance. We 
considered only those candidates that appear in at least three consecutive DMs to trace the 
DM dependence of signal-to-noise (SNR). A list of candidates was prepared after sorting. 


5.5. Candidate folding 

The zero-DM subtracted filter-bank data were dedispersed and folded with the period¬ 
icity parameters of sorted candidates in 32 sub-integrations and 64 frequency bands using 
PRESTO tool PREPFOLD. The left panel of Figure]^ shows an example of typical diagnostic 
plot from candidate folding for an MSP discovered in the GHRSS survey. 
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5.6. Candidate investigation 


A survey such as the GHRSS produces a large number of pulsar candidates, and it takes 
a lot of human effort to visually inspect the candidates in a reasonable amount of time. For 
example, a typical GHRSS survey scan with good RFI conditions produces around 100 to 
200 folded candidates which are then manually investigated for possible pulsed emission. 
However, in the case of a moderately bad RFI environment, the number of candidates can 
be a few thousand for each pointing. We employed a neural network based binary classiher 
capable of separating candidates arising from noise or RFI, from those generated by radio 


pulsars following Lyon et ah (2015). We used this scoring method and machine learning 


algorithm to hlter the best candidates for further inspection. The classiher is a very-fast- 
decision-tree and it uses the statistics (mean, standard deviation, kurtosis and skew) of the 
candidate prohle and how the signal-to-noise varies as a function of DM as input scores. The 
classiher is applied to all periodicity candidates to label those as positive or negative, and 
only the positive candidates are viewed by eye. For example, 7700 candidates from one epoch 
were reduced to just 156 positive predictions, resulting in about 50 times reduction in the 
number of candidates to be viewed. Out of these positive predictions, 17 were re-detections 
of known pulsars or their harmonics, and 5 were promising pulsar candidates for follow up. 
We also examined all of the same 7700 candidates by eye to conhrm that the classiher had 
performed accurately. The current version of the classiher was trained using a LOFAR data 
set but will be re-trained using the GHRSS data. The training set will include re-detections 
of known pulsars and discoveries from the GHRSS survey, as well as a sample of noise and 
RFI candidates from this survey. This will improve the performance of the classiher and 
reduce the number of false positives. 


5.7. Single pulse processing 


Besides the Fourier based periodicity searches, we performed a single pulse search to 
detect transient events, such as the RRATs or FRBs with possible extragalactic origin. We 
developed a PRESTO-based single pulse search pipeline for the Hydrus cluster of the Jodrell 
Bank pulsar group and IBM cluster at the NGRA for analysing the GHRSS survey data. 
After hlter-bank conversion, zero-DM subtraction and RFI bagging, we dedispersed the raw 
data hie at various trial DM values up to a DM of 2000 pc cm“^ following ddplan. In order 
to accelerate the single pulse search we reduced the nominal time resolution to 1 ms and the 
dedispersed time series is further down sampled according to ddplan. We then searched for 


single pulses in the time series using the presto code single.pulse.search.py (McLaughlin 


et al. 2006) above a detection signihcance of 5a. Then we investigated the diagnostic plots. 
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Figure [^shows output from our CPU based single pulse search pipeline for the known pulsar 
J1559—4438. For a detected dispersed event we further investigate the dynamic spectrum 
to confirm the expected quadratic time delay nature of an astrophysical signal. 


5.8. GPU-based processing 


Since analysis of our data were limited by available compute power, re-processing the 
data in the future as we upgrade our compute power will possibly result in additional dis¬ 
coveries. We are in the process of setting up a multi-node Graphics Processing Unit (GPU) 
cluster, consisting of 5-nodes each with 4 Nvidia GTX 980 GPUeach, at the University of 
Manchester to aid the GHRSS survey. Re-analysis of the GHRSS data in search for highly 
accelerated pulsars over a linear drift range of ± 250 frequency bins is being performed with 
a variant of the peasouf^^ software implemented by us for the GHRSS survey. In addition, 
a GPU based single pulse search pipeline based on heimdall ( Barsdell et al.||2012 ) software 
was optimised and implemented for the GHRSS survey. This integrated multi GPU pipeline 
for pulsar and transient search implemented for the GHRSS survey is called bieorst (paper 
in preparation). Figure]^ shows the output from this pipeline for single pulse searches. 


6. Survey sensitivity 


A pulsar is detectable (with 5a detection signihcance) in a survey with the incoherent 
array if it exceeds some minimum flux density (Fpuisar) that can be calculated using the 
radiometer equation: 


, ^ Tree + Tsky / W 

G^BN,N,t\l P-w 


( 1 ) 


where Tree and Tgky are the temperatures of the receiver and sky respectively, G the gain 
of individual antennas, B the bandwidth, A"p the number of polarizations, W the number 
of antennas, t the integration time, w the effective pulse width (including all instrumental 
smearing), and P the pulse period. As the telescope gain compensates for the larger band- 
widths available at high frequencies (~ Giow/Ghigh > ^/Bhigh/Blow), h makes more sense 
to conduct surveys at lower frequencies as the larger beam-width greatly increases the sky 
coverage (oc u~‘^) for a given dwell time (f) with similar sensitivity. 


^°http://www.geforce.co.uk/hardware/desktop-gpus/geforce-gtx-980/specifications 
^ ^ https: / /git hub .com/ewanbarr/peasoup 
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Pulsar period (ms) 



(a) 


(b) 


Fig. 7.— (a) Estimated sensitivity of the GHRSS (green), GBNCC (bine), HTRU (magenta) 
and GMRT-FERMI (cyan) snrveys. The solid lines are for DM 50 pc cm“^ and dashed lines 
for 100 pc cm“^. Sensitivities are scaled to 322 MHz with a spectral index of —1.7. (b) 
Estimated sensitivity of the GHRSS snrvey to MSPs as a fnnction of DM for a range of 
pnlse periods, 2 ms (green), 5 ms (bine), 10 ms (magenta) and 30 ms (cyan). 


For the GHRSS snrvey we calcnlate a theoretical survey sensitivity of ~ 0.5 mJy at 
322 MHz (considering the radiometer equation) for a 5a detection for 10% duty cycle, with 
the GMRT incoherent array of gain 2.5 K/Jy (~ G-^iVpiVa) for 32 MHz bandwidth and 


considering a system temperature of 106 K. Figure 7a plots the estimated sensitivity of 
the GHRSS survey compared with GBNGG at GBT and HTRU at Parkes and the Fermi- 
directed survey at GMRT at DMs of 50 and 100 pc cm“^. Figure T& illustrates that the 
GHRSS survey is the most sensitive survey in the GBNGG complementary sky. Figure 7^ 
plots the estimated sensitivity of the GHRSS survey for MSPs as a function of DM. This 
indicates that for a 5 ms pulsar even at a DM of 100 pc cm“^ the estimated sensitivity 
is ~ 0.75 mJy, which is only 1.5 times the predicted GHRSS sensitivity. Enhanced time 
and frequency resolution of GHRSS survey by a factor of 4 (Section]^ enables an increased 
search sensitivity for MSPs up to a factor of 3 compared to the Fermi-directed MSP survey 
at the GMRT, which used 512x0.0651 MHz channels sampled at 61.44 /rs. This is validated 
by observing J0248+4230 in these observing modes (Figure]^. With the newly implemented 
high resolution mode, the GHRSS survey is sensitive enough to detect MSPs at higher DMs 
as dispersion smearing is less than scatter broadening up to DM of 75 pc cm“^ (Figure]^. 
This will be useful considering 60% of the survey region will have DM< 100 pc cm“^ based 


on Gordes et ah (2002). Thus the GHRSS survey modes are optimal for the targeted survey 
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region. Figure 7a also illustrates that the GHRSS survey is the most sensitive survey in the 
GBNCC complementary sky. Table shows a comparison between the major off-Galactic 
plane surveys. 

A transient event of 5 ms duration will be detected at lOcr detection significance if it 
exceeds some minimum flux density (^transient) calculated from the radiometer equation: 




transient 


10 


T 

-L rf 


T. 


sky 


G^BNpN^T 


( 2 ) 


where r is the duration of the burst, Tree is 66 K and Tsky is 40 K. For the GHRSS survey 
we calculate a sensitivity of 1.6 Jy for a lOcx detection limit for 5 ms transient millisecond 


bursts, considering weak scattering (Thornton et ah 2013). 


7. Re-detection of known pulsars 


Re-detection of known pulsars and comparison with the catalogued flux density helps 
us to evaluate the survey sensitivity. The pulsars detected within the measurable primary 
beam gain pattern (i.e. up to the hrst null at ± 100') of the GMRT antennej^ were 
considered for flux density measurement. We had chosen only the highest signal-to-noise 
re-detection of a given pulsar, as some pulsars were re-detected in multiple adjacent scans. 
With 35% of the GHRSS survey complete, we had 23 such re-detections of known pulsars. 
There are 30 known pulsars above 15a detection significance present within the observed 
GHRSS survey pointings. Most of the non detections of known pulsars correspond to the 
pulsars being at the edge of beam and/or presence of RFI. We exploited the simultaneous 
availability of visibilities, to determine calibrated interferometric flux densities. Imaging was 
done using a pipeline that was developed specihcally for analysis of GMRT data (Kudale 
et ah, in preparation). The raw interferometric data were converted to FITS and then 


calibrated and flagged using the flagcal pipeline (Prasad et ah 2012; Ghengalur et ah 


2013). We performed flux calibrations using flux densities of the observed phase calibrators 
as given in the VLA calibrator manna The data were then imaged using a custom imaging 


pipeline built around utilities in the CASA package as well as the source identification package 
PYBDSM (paper in preparation), pybdsm was used to identify regions with emission in the 
dirty image which were set as the clean boxes in the cleaning and self-calibration cycles. 
First a wide, low resolution image was made (using mainly data from the central square 


^^http://www.ncra.tifr.res.in:8081/ngk/primarybeam/beam.html 
^^https://science.nrao.edu/facilities/vla/docs/manuals/cal 
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antennas). We nsed this to identify sonrces that lie ontside the half-power beam-width 
(HPBW) of the telescope. Then a high resolution image using all available antennas was 
made. The image covers the HPBW, but with extra “flanking fields” centred on all of the 
outlying sources identified in the low resolution image. The imaging and self calibration 
procedures typically included two rounds of phase-only self calibration followed by a final 
round of phase and amplitude calibration. Since most of the sources have been observed 
at very low elevation angles (as the fields in this survey were all southern fields) refractive 
effects can cause a signihcant shift in the position of the sources. Wherever possible this 
shift was calibrated using catalogued observations at higher frequencies (e.g. using the 843 
MHz SUMS^^ catalogue in helds that overlap with that survey, or sources with known 
positions in the NEDp^ database). After this correction we searched for emission around the 
known position of the pulsar. The measured flux density was then corrected for the primary 
beam attenuation to yield the final flux density of the pulsar. We could determine an 
interferometric flux density for 12 such re-detected fields having overlap with SUMSS/NED 
catalogued observations. Other re-detections could not be detected in the image due to 
increased noise near the pulsar position and/or significant flagging during imaging. Figure 
shows the re-detection of PSR J1312—5402 at 49' offset from the survey pointing centre, 
which clearly illustrates that pulsars of a few mJy can even be detected outside the HPBW 
of the survey pointing. Table lists the offset from the survey pointing centre, pulsed-SNR 
with primary beam correction and the corresponding imaging flux density of the re-detected 
pulsars. Similar to the GBNCC survey (Stovall et al. 2014) we estimated the expected 
catalogued flux density at 322 MHz for each of these pulsars, either by using the spectral 
index when available in the ATNF catalogue, or by extrapolating 400 MHz mean flux density 
reported in the ATNF catalogue assuming a spectral index of —1.7. Large uncertainties 
could be associated with the catalogued 400 MHz flux densities as well with the assumption 
of —1.7 spectral index, which indicates that the expected flux density estimates at 322 
MHz have low accuracy. Figure IR compares the observed interferometric flux density with 
the extrapolated catalogued flux density. The majority of the detections have observed 
flux densities within ±50% of the expected values (marked by green shaded region). In 
conjunction to the uncertainties involved in estimating the expected catalogued flux density, 
the observed sensitivity degradation for some of the re-detections could also be caused due 
to scintillation or significant flagging due to RFI. 


http://heasarc.gsfc.nasa.gov/W3Browse/radio-catalog/sumss.html 
^^https: //ned.ipac.caltech.edu 
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8. Discoveries 


In this paper we report pulsar and transient search results from 35% completion of the 
GHRSS survey covering 1000 deg^ of sky. With the periodicity search we discovered 10 
pulsars including one millisecond pulsar in a binary system. Figure [TT] shows the discovery 
profiles and Table lists the periods, dispersion measures and discovery flux density val¬ 
ues. Though some of the newly discovered pulsars are reasonably strong and might have 
been missed by previous surveys (e.g. J0418—4154, J0702—4956 and J0514—4407), we also 
discovered a few weak pulsars, indicating the highly sensitive nature of the survey even in 
the presence of some residual RFI after flagging. At the discovery epoch the newly dis¬ 
covered MSP had an estimated flux density of 1.6 mJy which is close to the theoretical 
sensitivity considering the system temperature of 109 K (Haslam et ah 1982) and 30% duty 
cycle. We confirmed all the discoveries in subsequent epochs of observations. The GHRSS 
discoveries have dispersion measures ranging from 15.4 pc cm“^ (for PSR J0514—4407) up 
to 133 pc cm“^ (for PSR J1456—48) and periods from 5.04 ms (for PSR J2144—5237) up to 
1169.89 ms (for PSR J1559—44). Most of the newly discovered pulsars (e.g. J0418—4154, 
J0702-4956, J0919-42, J1255-46, J1456-48, J1559-44 and J1708-52) have single peaked 
profiles. PSR J1947—43 has a relatively wide profile and there is a hint of a double peaked 
profile. For PSR J0514—4407 a double component pulse profile is observed. The profile 
shapes were confirmed with multi-epoch observations. The 52.04 ms rotational period of 
PSR J1255—46 puts it in the domain of young pulsars or double neutron star relativistic bi¬ 
naries. A more detailed study of this pulsar is in preparation. The 5.04 ms MSP J2144—5237 
has a multi-component prohle. The right panel of Figurej^shows the sensitive coherent array 
observations aided by precise localisation (Section]^. The flux density of this MSP varies 
dramatically due to scintillation, indicating the need for re-observations of candidate MSPs 
with longer integrations. A drift over Fourier frequency by 4x10“^ Hz, which corresponds 


to a period derivative of —1.4x 10“^^ s s“^ or a line-of-sight acceleration of 1.4 m s“^ ob¬ 
served for J2144—5237 over two hours, indicates that this MSP is in binary system. PSR 
J2144—5237 could be a good timer for inclusion in the International Pulsar Timing Array. 


The area normalised discovery rate of the GHRSS survey is 0.01 pulsars per deg^ which 
is higher than the 0.006 pulsars per deg^ achieved with the GBNGG survey (Section |^. 
It may be that this higher discovery rate is purely due to low number statistics at this 
stage or could perhaps represent an underestimate of the sensitivity of the GMRT in the 
incoherent mode or perhaps an overestimate of the GBNGG sensitivity. It may be useful to 
survey an overlapping region of sky to enable a more direct comparison between these low 
frequency surveys. However the GHRSS pulsar discovery rate is lower than the prediction 
by PSRPOPPY prediction (Section like many of the other surveys. 
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We have not detected any FRB or RRAT with the single pulse processing in the GHRSS 
survey. However, considering the FRB rate predicted in Section we should have detected 
only about 2~^l FRBs in this survey region based on a flat spectral index consideration. 
Absence of which is not constraining enough considering factor of few lower FRB prediction 
by recent studies (Rane et al. 2015 Champion et al. 2015; Petroff et ah 2014). 


9. Localisation of the GHRSS discoveries 

We utilised the interferometric imaging capability of the GMRT to localise pulsars and 
transient events. Such localisation in the image plane is very important for FRBs for host 
galaxy identification and multiwavelength follow up. 


9.1. Localisation of the pnlsars with gated correlator 


A gating interferometer is an excellent tool to improve the noise statistics in order to 
increase the detection significance towards a pulsed signal. The development of a coherently 


dedispersed gated correlator for the GMRT is detailed in Roy et al. (2013). In this design, 


the visibility time-series were derived off-line from the recorded raw voltage data and were 
binned using multiple gates and on each gate they were folded using the best-fit topocentric 
rotational model, derived from the synchronously generated incoherent beam data. We 
improved the noise statistics further by subtracting the visibilities of adjacent gates. The 
ON-OFF subtracted visibilities were then calibrated and flagged using flagcal pipeline 


(Prasad et al. 2012 Ghengalur et al. 2013). The ON-OFF image was made using the AlPS 
package, where the pulsar is unambiguously localised as a point source. For pulsars with 
10% duty cycle, we expect to get a factor of 3 improvement in detection significance in the 
gated image plane compared to the continuum plane. The positional uncertainty therefore 
reduces from 80' (half-power beam-width at 322 MHz) to <20" (typical synthesized beam 
used in the image made at 322 MHz). 

We recently started follow up of these pulsars using the coherently dedispersed gated 
correlator which resulted in the successful localisation of 4 pulsars, namely: J0418—4154, 
J0702—4956, J0514—4407 and J2144—5237 with an accuracy of <10" (positions are listed 
in Table . Figure plots the ON-OFF gated images of these 4 pulsars. Such timing 
independent, rapidly determined and fairly precise positions allow us to use sensitive coherent 
beams (up to 5x the sensitivity of the incoherent beam) of the GMRT for follow up timing 
























23 


observations, substantially reducing the use of telescope time. Further improvement of the 
detection signihcance is seen while using the coherent array, as it is better resistant to RFI. 
Right panel of Figure shows the huge improvement in sensitivity achieved with coherent 
array observations for PSR J2144—5237. Moreover, such a priori astrometric models allow. 


rapid convergence in pulsar timing as discussed in Section 10 


9.2. Localisation of FRBs in the image plane 


In addition to localisation of the pulsars, simultaneous availability of the visibility data 
during the high resolution survey provides us the opportunity to localise the transient events 
detected in time-domain searches. All 17 known FRBs have been discovered with single 
dish telescopes and interferometric detection of an FRB is of signihcant importance. This 
technique is successfully utilised in the study of the eclipsing pulsars with the GMRT (e.g. 
PSR J1227—4853 (Roy et ah 2015)). An example of interferometric localisation of a re¬ 


detection of a known pulsar in the GHRSS survey (Figure also illustrates the pipeline 
that will be followed for localisation of possible transient detection in the GHRSS survey. 
Assuming that interferometric gain is similar to coherent array gain (i.e. 5 times of incoherent 
array), FRBs with flux density > 3 Jy could be localised at > 5a signihcance in the image 
plane from simultaneously recorded visibilities at 2 s resolution. 


10. Timing and follow up of the GHRSS discoveries 


We are continuing regular timing follow up of the GHRSS pulsars, with the coherent 
array of the GMRT, to determine their rotational and binary parameters (for the MSP) and 
place them on the P-P diagram. Long term timing will also allow us to probe the possible 
optical and y-ray counterparts of the newly discovered pulsars. The timing-independent 
rapid astrometric measurement of newly discovered pulsars using gated imaging helps to 
break the degeneracy in the timing £t caused by co-variance between pulsar position and 
rotational parameters. We successfully demonstrated this by timing two of the localised 
pulsars, J0418—4154 and J0702—4956, using sparsely spaced TOAs over a span of ~ 200 
days. Gonsidering reasonable fractional bandwidth of 10% at 322 MHz, we determined the 
DM using timing fits of sub-band TOAs. The timing ephemeris of both the pulsars are given 
in Table[^ The post-£t residuals are shown in Figurejl^ Residuals for J0702—4956 are fitted 
with three timing models having different astrometric positions: the gated interferometric 
position, 1' offset and 5' offset from the gated interferometric position, which is illustrated 


in Figure 14 The residuals increase as the position offset increases and the residuals are 
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wiggled in phase at 5' offset. This warrants more closely spaced and regularly sampled 
TOAs to overcome the phase ambiguity caused by position error. Thus a priori knowledge of 
a pulsar position from the gated-image is extremely important, considering that the discovery 
position can be anywhere within the HPBW (± 40') or even outside the HPBW (e.g. Figure 

0 . 

Following the localisation of a GHRSS pulsar we routinely check for nearby sources 
from existing opticap^ and high-energy catalogue^^ For one of the GHRSS discovery pul¬ 
sar, J0514—4407, we found a possible counterpart 3FGL J0514.6—4406, with a position 
uncertainty of about 0.16°. This is only 1.8' from the pulsar, so well within predicted y-ray 
position. However, none of the other pulsars listed in Tablej^have 2FGL/3FGL counterparts. 
We plan to fold the Fermi photons with the timing model of the newly discovered pulsars to 
check if any of these pulsars exhibit y-ray pulsations. This will provide input towards the 
radio-loud versus radio-quiet y-ray pulsar ratio, which is a reasonable differentiator between 
y-ray population models. 


11. Relevance to the SKA 


Discoveries of pulsars, in particular MSPs, are hindered by their radio faintness that 
demands deeper searches with larger telescopes and we are reaching the limit of what is 
possible for fully steerable single dishes. Large arrays of many smaller telescopes are the 
future for large radio telescopes leading ultimately to the world’s largest telescope the SKA. 
The SKA’s wide field-of-view, high sensitivity, multi-beaming and sub-arraying capabili¬ 
ties, coupled with advanced pulsar search backends, will result in the discovery of a large 


population of pulsars (Keane et al. 2015). As it is currently the largest telescope in the 


metre-wavelength regime, the GMRT, an SKA pathhnder telescope, is the prototype for the 
SKA in many ways, and provides an excellent test bed for new techniques. Tools for efficient 
RFI mitigation, optimised search techniques, synchronous beam-formation and imaging will 
provide vital input to the SKA. In particular the high resolution GHRSS data are being 
used for testing the optimised dedispersion and de-acceleration plan of the pulsar search 
sub-element (PSS) of SKA. The periodicity search for the GHRSS needs an 8 million point 
transform, similar to the requirement of the SKA. Moreover, the rapid convergence in pulsar 
timing model possible with precise localisation of the pulsars with gated imaging described 
in Section]^ and IR demonstrates a feature that may be applicable in the SKA. 


ht tp: / / simbad. u—st rasbg. fr / simbad / sim—fcoo 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/4yr .catalog/ 
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The GHRSS survey is at a frequency that overlaps with the SKAl low frequency instru¬ 
ment (SKAl—Low) and considering SKAl-Low will have sensitivity ~ 0.05 mJy at 350 MHz 
(Keane et al. 2015) (i.e. 10 times GHRSS sensitivity) all the pulsars discovered with the 
GHRSS survey will be able to be timed and studied with SKAl—Low, providing an excellent 
opportunity to prepare for searches and timing for this next generation telescope. The GPU 
based processing pipeline implemented by us, biforst, performs a time-domain acceleration 
search on pulsar survey data, and the processing steps roughly follow those of the PSS time- 
domain processing pipeline. The GHRSS survey data are concurrently being processed by 
time-domain acceleration search using GPU and frequency-domain acceleration search using 
GPU, and thus provides detailed comparison between the algorithms, which can be a vital 
input to SKA choices to decide the appropriate search methodology. These results will pro¬ 
vide very valuable input into directing the prototyping and design effort for pre-construction 
stage. In future for the GHRSS survey we plan to use the baseband recording facility of 
the GMRT Software Backend followed by multi-pixel coherent beamformation (Roy et ah 
2012). Pixelisation of the full £eld-of-view while scanning the extreme southern sky with 
the GMRT will provide 3x sensitivity improvement and 20x increase in compute cost. For 
the extreme southern sky we will require 300 beams which is ~ 10% compute cost of SKA 
pulsar search. 


12. Conclusion 

In this paper we described the system configuration and initial discoveries from the 
GHRSS survey. We described the high resolution observing modes developed for this survey 
and the data processing pipeline for periodicity and single pulse search. Our survey sensitiv¬ 
ity is comparable to other ongoing low frequency surveys, and we are targeting a portion of 
sky which has not been surveyed at lower frequencies like 322 MHz for the last two decades. 
With 35% completion of the GHRSS survey covering 1000 deg^ of sky we discovered 10 
pulsars including 1 MSP which corresponds to one of the best pulsar per deg^ discovery 
rates of any survey off the Galactic plane. We also re-detected 23 previously known pulsars 
with the expected detection significance, which were in-beam pulsars, not necessarily at the 
pointing centre. The simultaneous time-domain search and imaging capability of this survey 
provides an opportunity of discovery and localisation of pulsars and transients. Utilising this 
we localised 4 of the newly discovered pulsars in the gated image plane. With the aid of ~ 
10" localisation, we obtained the timing solutions of 2 of the discoveries with relatively short 
data spans (< 1/2 year). In addition to the regular emission from pulsars, the GHRSS survey 
will reveal transient events like the RRATs or the FRBs. With an ongoing enhancement of 
processing power we expect more results soon. 
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Table 1: Survey parameters of the GHRSS survey for the MGL and HGL components 


Survey 

MGL 

with HiResl 

HGL 

with HiRes2 

Galactic region 

5<|6|<20° 

|6| > 20° 

Declination 

-40°<Dec<-54° 

-40°<Dec<-54° 

Integration time 

1200 s 

900 s 

Sampling time 

61.44 fis 

30.72 /rs 

Bandwidth 

32 MHz 

32 MHz 

Number of channels 

2048 

1024 

Frequency Resolution 

15.625 kHz 

31.25 kHz 

Number of pointings 

682 

911 

Sky coverage 

1227 deg2 

1639 deg2 

Data/pointing 

37 GB 

28 GB 

Total data 

25 TB 

25 TB 
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Table 2: Major ongoing or recently completed off-Galactic plane surveys 


Survey name 
— Telescope 

Frequency 
of search 
(MHz) 

Sky coverage 

Discoveries 

Sensitivity! 

(mJy) 

HTRU^ 

1352 

-120°<1, 1<30° 

104 PSR, 26 MSP 

1.5 



|6|<15° 



— Parkes 


4500 deg^ 



HTRU-N 

1360 

l&l > 15°, Dec > -20° 

12 PSR 

1.5 

— Effelsberg 





GBNGG^ 

350 

Dec>—40° 

108 PSR, 12 MSP 

0.6 

- GBT 


19500 deg2 



GBTdriftscan^ 

350 

-21° <Dec <26° 

26 PSR, 7 MSP 

0.9 

- GBT 





A0327^ 

327 

0° <Dec< 28° 

24 PSR, 3 MSP 

0.3 

— Arecibo 





LOTAAS^ 

135 

Dec >0 ° 

18 PSR 

0.3 

-LOFAR 


5156 deg2 



GHRSS^^^ 

322 

-20° <Dec < -54° 

10 PSR, 1 MSP 

0.5 

- GMRT 


1000 deg^ 




f: Sensitivity calculated at 322 MHz considering 5a limit for 10% duty cycle for spectral 
index of —1.7 

Full GHRSS survey will cover 2866 deg^. 

http://researchbank.Swinburne.edu.au/vital/access/manager/Repository/swin;31985 

http://arcc.phys.utb.edu/gbncc/ 

http://astro.phys.wvu.edu/GBTdrift350/ 

http://www.naic.edu/deneva/drift-search/ 

http://www.astron.nl/lotaas 

http://www.jb.man.ac.uk/research/pulsar/Resources/ghrss.html 
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Table 3: Parameters for psrpoppy simulations 


Parameters 

Normal pulsars 

Millisecond pulsars 

Spin Period Distribution 

Log-normal 

Lorimer et al.j (2012)1 

logiQ{V (ms)) 

2.7 


std logiQ^P (ms)) 

0.34 


Luminosity Distribution 

Log-normal 

Power-law 


{logio L (mJy kpc^ )) = -1.1 

L(mJy kpc‘^)min = 0.1 


std(/o5fio L (mJy kpc^ ))= 0.9 

L(mJy kpc^)^^^; = 1000 
pow=—1.45 

Galactic Scale Height 

Exponential 

Exponential 


0.33 kpc 

0.5 kpc (Levin et al. 2013) 

Radial Distribution model 

Lorimer et al. (2006) 

Lorimer et al. (2006) 

Galactic electron distribution 

NE2001 model 

NE2001 model 


(Cordes et al. 2002) 

(Cordes et al. 2002) 

Spectral index distribution 

Gaussian 

Gaussian 

(«) 

-1.7 

-1.7 

std (a) 

0.35 

0.35 


f: with this built-in option in PSRPOPOY, periods of millisecond pulsars are picked in at 
random from data base provided in Lorimer et al. (2012). 


Table 4: Dispersion plan and corresponding time resolution for HGL and MGL according to 
the DDPLAN of PRESTO 


DM range 

DM step (pc cm ^) 

HGL 

Sampling resolution (ps) 

DM step 

MGL 

Sampling resolution (ps) 

1-25 

0.010 

30.72 

0.010 

61.44 

25-50 

0.025 

61.44 

0.010 

61.44 

50-100 

0.050 

122.88 

0.025 

61.44 

100-150 

0.100 

245.76 

0.050 

122.88 

150-300 

0.200 

491.52 

0.100 

245.76 

300-500 

0.500 

491.52 

0.200 

491.52 
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2 Pulses of Best Profile 



Candidate: PSR_024a+4230 

Telescope: GMRT 

Epochj^ = 56942.74-359360817 

Epochj, - N/A 

Tsamp.» = 6.144e-05 

Data Folded - 29368320 

Dota Avg = 5.089e-f04 

Dato StdDev = 228.6 

Profile Bins = 64 

Profile Avg = 2.335e+10 

Profile StdDev = 1.549e+05 


Search Information 

RAjjq3P = 02:48:29.8900 DEC^j^qoo = 42:30:13.9400 

Folding Parometera 

D0F,p = 32.27 = 3.304 P(Noise) < Z.06e-17 (8.4ff) 

Dispersion Measure (DM: pc/cm^) — 48.250 

(ms) = N/A 


Ptcpo = 2.60070577(25) 
P'lopo (s/a) - 0.0(1.1)x10-’* 
'^"lopo (®/s ) = 0.0(3.9)x10' 
Binary Parometera 
Por# (=) = N/A 
OlSln(i)/c (s) = N/A 

- N/* 


^bory • 

P'bar, (s/s) - N/A 
P-'j; (s/a^) = N/A 

e = N/A 
u (rad) = N/A 



2 Pulses of Best Profile 



Candidate: PSR_0248-l-4230 
Telescope; GMRT 
Epochj^ = 56942.76958570489 
^P®'=^b«rv ■ N/A 


Dota Avg 
Data StdDev = 110.9 

Profile Bins = 64 

Profile Avg = 1.Ole- 

Profile StdDev = 7.53e- 


6.144e-05 

29490688 

2.191e+04 


Search Information 
3 = 02:48:29.8906 DECj^, 

Folding Porometere 


= 42:30:13.9375 


DOF,^ = 32.27 = 6.905 P(Noige) < 5.33e-57 Cl5.9a) 

Dispersion Measure (DM; pc/cm^) — 48.250 

• • (ms) » N/A 


3 (ms) = 2.60070662(13) 
('/a) - 0.t)(5.6)«10-'' 
P"i.,. (s/a ) = 0.0(2.0)i.l0- 
Binary Parameters 
P« (a) = N/A 
o,sin(r)/c (s) = N/A 
Tpea = N/A 


P'bar, (s/s) - N/A 

P-bJi (s/s") = N/A 

e = N/A 
w (rod) - N/A 



Phose 


Reduced ^ 


DM (pc/cm=) 


Period - 2.60070862 (ms) 


Fig. 8. — Output of PRESTO based GHRSS search pipeline for the MSP J0248+4230 to 
demonstrate the improvement in detection significance with the high resolution observing 
modes developed in the GHRSS survey. Top panel: with 32 MHz observing band split into 
512 frequency channels recorded every 61.44 /is, Bottom panel: with 32 MHz observing band 
split into 2048 frequency channels recorded every 61.44 /xs. 






























34 



ighj^m ^Qm Q^m 

Right Ascension (J2000) 


Fig. 9.— Pulsar J1312—5402 was detected outside the HPBW (marked by bigger circle) 
and 49' offset (marked by smaller circle) from the survey pointing centre with 17a pulsed 
detection significance. The interferometric flux density is 22 mJy. The right plot shows a 
small portion of the image containing the pulsar. 


Table 5: Parameters of the pulsars discovered in GHRSS survey 


Pulsar name 

Period 

(ms) 

Dispersion measure 
(pc cm“^) 

Detection significance 

Flux density! 
(mJy) 

PSR J0418-4154 

757.11 

24.5 

50 

10.3 

PSR J0514-4407 

302.2 

15.4 

42 

9.7 

PSR J0702-4956 

666.66 

98.7 

30 

15.7 

PSR J0919-42 

812.6 

57.8 

19 

6.4 

PSR J1255-46 

52.0 

42.9 

12 

0.8 

PSR J1456-48 

536.81 

133.0 

15 

1.2 

PSR J1559-44 

1169.89 

122.0 

8 

1.7 

PSR J1708-52 

449.62 

102.6 

9 

1.4 

PSR J1947-43 

180.94 

29.9 

17 

4.7 

PSR J2144-5237 

5.04 

19.0 

9 

1.6 


1 : Flux density is without primary beam correction. 


iJy/Beam 
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Fig. 10.— Comparison of catalogued mean flux density at 322 MHz (extrapolated either by 
using the spectral index when available in the ATNF catalogue, or extracted considering 400 
MHz mean flux density reported in the ATNF catalogue and a spectral index of —1.7) and 
observed interferometric flux density at 322 MHz for 12 re-detected pulsars. The shaded part 
(green) represents a region with the observed flux density within ±50% of the catalogued 
flux densities. 
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(f) J1456-48 


(g) J1559-44 


(i) J1947-43 (j) J2144-5237 


Fig. 11.— Discovery plots of the ten pulsars discovered in the GHRSS survey. The observa¬ 
tions were performed in the incoherent array mode. The folded prohles are plotted (twice) 
in the top panel of each plot. Signal strength as a function of rotational phase is plotted in 
the bottom panel. Severe RFI occasions are masked, however as can be seen in some cases 
signihcant RFI remains. 





Table 6: List of 23 re-detected pulsars in the GHRSS survey. Offsets from the GHRSS 
beam centre, primary beam corrected pulsed SNR and interferometric flux densities (for 12 
pulsars, the others either did not had overlapping SUMSS/NED observations, or could not 
be detected in the image due increased RMS near the pulsar position and/or signihcant 
flagging during imaging) are tabulated. 


Pulsar name 

offset (') 

Pulsed-SNR 

Interferometric Flux density (mJy) 

J0255-5304 

25 

31 

- 

J0745-5353 

72 

272 

- 

J0809-4753 

69 

176 

- 

J0842-4851 

53 

60 

- 

J0840-5332 

53 

30 

- 

J0843-5022 

0.0 

11 

12±1 

J0905-5127 

9 

28 

- 

J0907-5157 

51 

130 

72±3 

J0924-5302 

26 

77 

26±1 

J0934-5249 

54 

46 

60±10 

J0955-5304 

62 

132 

- 

J1003-4747 

2.5 

82 

11.0T0.9 

J1036-4926 

43 

52 

6.0±2 

J1045-4509 

33 

50 

24±1 

J1143-5158 

40 

64 

- 

J1240-4124 

36 

93 

23±2 

J1312-5402 

58 

43 

41±1 

J1320-5359 

79 

428 

20±2 

J1355-5153 

47 

114 

50±5 

J1544-5308 

15 

36 

- 

J1559-4438 

3.7 

59 

120±5 

J1902-5105 

75 

159 

- 

J2241-5236 

13 

71 

- 
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Fig. 12.— Localisation images for 4 pulsars discovered in the GHRSS survey: J0418—4154, 
J0514—4407, J0702—4956 and J2144—5237. The offsets of the pulsar from the pointing 
centre is given in Table Since J0514—4407 and J2144—5237 are at large offsets from the 
pointing centre, we show a smaller held closer to the pulsar. 







































Table 7: Parameters obtained from the gated imaging 


PSR 

Gated J2000 
position 
(") 

Offset from 
pointing 

centre 

Nnmber 
of gates 

Observing 

dnration 

(s) 

Gated 
flnx density 
(mJy) 

Gated 

SNR 

J0418-4154 

04‘'18“02!88(T.'4); 

-41°54'1T.'89(7"8) 

7.4' 

24 

900 

88 

9 

J0514-4407 

05*"14“5P84(T.'04); 

-44°07'06'.'51(8'.'4) 

32' 

21 

900 

20 

5 

J0702-4956 

07^'02“38!54(T.'2); 

-49°56'56'.'27(7"8) 

6.2' 

16 

900 

29 

6 

J2144-5237 

21M4“39!2(0'.'7); 

-52°37'32'.'17(3'.'8) 

28' 

10 

3600 

6 

12 


fO 




322 MHz TOAs ' 


100 150 200 250 

MJD - 56900 


(a) 


6 - 


0 - 


-6 - 


£-10 


(b) 


322 MHz TOAs ' 


80 100 120 140 

MJD - 56900 


Fig. 13.— Post-fit timing residnals for J0418—4154 (left panel) and J0702—4956 (right 
panel). Mnltiple TOAs are derived for each observing epoch. Post-£t residnal is 1.1 ms for 
J0418-4154 and 1.0 ms for J0702-4956. 
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Fig. 14.— The post-fit timing residual for PSR J0702—4956 fitted with three timing models 
having different astrometric positions, the gated interferometric position (solid red line), 
V offset (green dashed line) and 5' offset (blue dotted line) from the gated interferometric 
position. The errors on TOAs are magnified by 10 x for better visualisation. A ±0.5 phase 
wrap is noticed for the timing residual with 5' position offset. 










Table 8: Timing parameters of PSR J0418—4154 and J0702—4956 


Parameters 

J0418-4154 

J0702-4956 

Gated imaging position 

Right ascension(J2000). 

04^18“02!88±R4 

07^02“38!54±R2 

Declination(J2000). 

-41°54'H"89±A'8 

-49°56'56"27±A'8 

Parameters from radio-timing" 

Right ascension(J2000). 

04^18“04!2(2) 

07‘"02“39!3(2) 

Declination(J2000). 

-41°54'10'.'8(6) 

-49°56'34"2(6) 

Pulsar frequency /(Hz). 

1.3207965590(1) 

1.501515953(1) 

Pulsar frequency derivative / (Hz s“^) 

-1.7 (l)xlO-i® 

0.00 

Period epoch(MAD). 

57055 

57086 

Dispersion measure DM(pc cm“^).... 

24.54(8) 

98.7(1) 

Span of timing data(MAD). 

56915.88-57180.25 

56942.95-57086.64 

Number of TOAs. 

44 

29 

Post-fit residual rms(ms). 

1.1 

1.00 

Reduced chi-square. 

1.05 

0.7 


Errors in the last digit are in parentheses 














